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We observe evidence for CP violation in the decay D"*" Ks'"''^ using a data sample with an 
integrated luminosity of 977 fb~^ collected by the Belle detector at the KEKB e^e~ asymmetric- 
energy collider. The CP asymmetry in the decay is measured to be (—0.363 ± 0.094 ± 0.067)%, 
which is 3.2 standard deviations away from zero, and is consistent with the expected CP violation 
due to the neutral kaon in the final state. 
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In the standard model (SM), violation of the combined 
charge-conjugation and parity symmetries (CP) arises 
from a nonvanishing irreducible phase in the Cabibbo- 
Kobayashi-Maskawa flavor- mixing matrix [ij. In the SM, 
CP violation in the charm sector is expected to be very 
small, 0(0.1%) or below 0. Since the discovery of the 
J/ip 0] and the subsequent discovery of open charm 
particles [3], CP violation in charmed particle decays 
has been searched for extensively and only recently be- 
came experimentally accessible. To date, after the FO- 
CUS [5], CLEO [6], Belle 0, and BaBar @ measure- 
ments, the world average of the CP asymmetry in the de- 
cay D+ KItt+ fgl] is (-0.54±0.14)%, which is the first 
evidence oi CP violation in charmed particles. However, 
it should be noted that the observed asymmetry is con- 
sistent with that expected due to the neutral kaon in the 
final state and is not ascribed to the charm sector. Re- 
cently, LHCb reported AAcp = (-0.82 ± 0.21 ± 0.11)%, 
where /S.Acp is the CP asymmetry difference between 
D° K+K- and I?" tt+tt- decays This is 

the first evidence of non-zero AAc p in charmed particle 
decays from a single experiment. 

In this Letter we report the first evidence for CP vio- 



lation in charmed meson decays from a single experiment 
and in a single decay mode, D'^ — )■ K^tt^ , where de- 
cays to TT+TT". The CP asymmetry in the decay, Aqp, 
is defined as 
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^CP 



the partial decay width, and A'^p and 
denote CP asymmetries in the charm decay 



mixing in the SM [12, ll3[ , respec 



(AC) and in K° 
tively. 

of amplitudes for — )■ K^Tr~^ and — )■ K'^tt^ decays 
where the former is Cabibbo-favored (CF) and the latter 
is doubly Cabibbo-suppressed (DCS). In the absence of 
direct CP violation in CF and DCS decays (as expected 
within the SM), the CP asymmetry in D 
decay within the SM is A\ 

(-0.332±0.006)% [11 from semileptonic decays [l5|. 
On the other hand, if processes beyond the SM con- 
tain additional weak phases other than the one in the 
Kobayashi-Maskawa ansatz [H] , interference between CF 



The observed ^^^571'"'" final state is a coherent sum 



^p, which is measured to be 
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and DCS decays could generate an 0(1%) direct CP 
asymmetry in the decay — )■ K^tt'^ [l^. Thus, ob- 
servation of Ac p inconsistent with A^^p in D+ — >■ K^-k'^ 
decay would be strong evidence for processes involving 
new physics [13l Il6|. 



We determine A(~,p 
try in the signal yield 



by measuring the asymme- 



a; 



(2) 



where A^rcc is the number of reconstructed decays. The 
asymmetry in Eq. ([2]) includes the forward-backward 
asymmetry (Afb) due to j*-Z^ interference and higher 
order QED effects in 



e cc [i/IJ, and the detec- 
tion efficiency asymmetry between tt"*" and 7r~ {A'^^) as 
well as Acp- In addition, Ref. [l^ calculated another 
source denoted due to the differences in interactions 
of K'^ and mesons with the material of the detector. 
(The existence of this effect was pointed out in Ref. [7].) 
Since we reconstruct the Kg with tt+tt^ combinations, 
the 7r+7r~ detection asymmetry cancels out for Kg. The 
asymmetry of Eq. © can be written as 



Ai 



+ (plpt+^cos^j?) + Ar,(p'^\) (3) 



by neglecting the terms involving the product of asym- 
metries. In Eq. ([3]), Acp is independent of all kinematic 
variables other than K^ decay time due to the K^ in 
the final state [31 , ^fb odd function of the cosine 
of the polar angle of the Z?"*" momentum in the center- 
of-mass system (CMS), A^ depends on the transverse 
momentum and the polar angle of the 7r+ in the lab- 
oratory frame (lab), and Ax> is a function of the mo- 
mentum of the K'g in the lab. To correct for Af^ in 
Eq. ©, we use D+ K-ti+tt+ and D° K-jt+tt^ 
decays, and assume the same ApB for Z?+ and 
mesons. Since these are CF decays for which the di- 
rect CP asymmetry is expected to be negligible, in anal- 



ogy to Eq. ©, A 



D+^K-TT+TT + 



and At, 



D" 



include 

+ 



one can measure A^* by sub- 



ApB, Af , and . Thus with the additional 
term in A^^^^ 

tracting A^"^-^' ' ^" from A'^^'^^'^ We obtain 

Ad according to Ref. [18j]. Using ^roc s'^corr sJ^q^^^j^ 

in Eq. Q, which is "^^^'^ after the A^^ and Ad 
corrections. 
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Afsicose^T), (4) 



we extract Acp and Afb using 



A 



CP 



= [Ai 

+ Ai 



■^FB 



"(+COS0g^S) 

"(-cos0™S)]/2, (5a) 
"(+cos0g^S) 



- Afe:^''-^'^-"(-cos0g^«)]/2. (5b) 



Note that extracting Acp in Eq. ([4]) is crucial in Belle 
due to the asymmetric detector acceptance in cos 6*^^^ . 

The data used in this analysis were recorded at the 
T{nS) resonances (n = 1,2,3,4,5) or near the T(45') 
resonance with the Belle detector [23| at the e^e~ 
asymmetric-energy collider KEKB [2l| . The data sample 
corresponds to an integrated luminosity of 977 fb~^. 

We apply the same charged track selection criteria that 
were used in Ref. without requiring associated hits 
in the silicon vertex detector [23|. We use the stan- 
dard Belle charged kaon and pion identification . We 
form Kg candidates from tt+tt" pairs, fitted to a com- 
mon vertex and requiring the invariant mass of the pair 
Min+TT-) to be within [0.4826,0.5126] GeV/c^, regard- 
less of whether the candidate satisfies the standard Kg 
requirements jS^]. (We refer to the Kg candidates not 
satisfying the standard criteria as "loose Kg" .) The Kg 
and TT"*" candidates are combined to form a candi- 
date by fitting them to a common vertex and the Z?^ 
candidate is fitted to the e"'"e~ interaction point to give 
the production vertex. To remove combinatorial back- 
ground as well as mesons, which are produced in 
possibly CP violating B meson decays, we require the 
Z?+ meson momentum calculated in the CMS (p^+) to 
be greater than 2.5 and 3.0 GeV/c for the data taken 
at the T(4S') and T(5S') resonances, respectively. For 
the data taken below T(4S'), which is free of B mesons, 
we apply the requirement p^+>2.0 GeV/c. In addition 
to the selections described above, we further optimize 
the signal sensitivity with four variables: the of the 
Z3+ decay and production vertex fits (x|> and x%), the 
transverse momentum of the tt"*" (pT7r+)i ^Lud the an- 
gle between the momentum vector, as reconstructed 
from the daughters, and the vector joining the -D"*" pro- 
duction and decay vertices (^) An optimization is 
performed by maximizing Af.q/y/JVfi + Mb with the four 
variables varied simultaneously [25| , where Afs + Mb and 
Mb are the yields in the Ks^~^ invariant mass signal 
([1.855,1.885] GeV/c^) and sideband ([1.825,1.840] and 
[1.900, 1.915] GeV/c^) regions, respectively. The optimal 
set of (xz)j Xpj Pttt+j requirements are found to be 
(<100, <10, >0.50 GeV/c, <160°), (<100, <10, >0.45 
GeV/c, <170°), and (<100, <10, >0.40 GeV/c, no re- 
quirement) for the data taken below the T(45'), at the 
T(4S'), and at the T(5S'), respectively. The candi- 
dates with the loose Kg requirement are further opti- 
mized with two additional variables which are the of 
the fit of pions from the Kg decay and the pion from the 
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D+ meson decay to a single vertex (x3,r)' ^^'^ the an- 
gle between the Kg momentum vector, as reconstructed 
from the daughters, and the vector joining the and 
Kg decay vertices (C). The two variables are again varied 
simultaneously and the optimum is found to be X3^>6 
and C<4° for all data. The inclusion of candidates 
with the loose Kg requirement improves the statistical 
sensitivity by approximately 5%. After the final selec- 
tions described above, there remains a background with a 
broad peaking structure in the Kgir'^ invariant mass sig- 
nal region, due to misidentification of charged kaons from 
Df — )• K^K~^ decays. The _D+ 7r+7r~7r+ background 
is found to be negligible from simulation [26| . Figure [T] 
shows the distributions of M{K'^Tr+) and M{K'^7r-) to- 
gether with the results of the fits described below. 




1.8 „ 1.9 



M(K^7t:^) (GeV/c-^) 



ccr^lO° 
o 



CD 



to 



CD 



10^ 




M(K^7t:-) (GeV/cO 



FIG. 1: Distributions of M(A'g7r+) (top) and M{K'^tt-) (bot- 
tom). Dots with error bars are the data while the histograms 
show the results of the parameterizations of the data. Open 
histograms represent the Z?^ Kgn^ signal. Shaded and 
hatched regions are KgK^ misidentification and com- 

binatorial backgrounds, respectively. 



neous fit to the D+ and D~ candidate distributions. Be- 
sides the asymmetry and the total signal yield, the com- 
mon parameters in the simultaneous fit are the widths of 
the Gaussian and the bifurcated Gaussian and the ratio 
of their amplitudes. The asymmetry and the sum of the 
D+ and D" yields from the fit are (-0.146±0.094)% and 
(1738±2)x 10'^, respectively, where the errors are statis- 
tical. 
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FIG. 2: map in bins of pi^^ and cos S'^'^ of the tv+ ob- 

tained with the D'^ K~tv'^tv'^ and — > K~tt^tv'^ samples 
(triangles). The A^^^ ^ map is also shown (rectan- 

gles). 



The K%-K^ signals are parameterized as a sum 

of a Gaussian and a bifurcated Gaussian distribution 
with a common mean. The combinatorial background is 
parameterized with the form e"+/5*f(^s'^ )^ where a and 
/3 are free parameters. The shapes and normalizations of 
the — > K^gK^ misidentification backgrounds are ob- 
tained with taking the asymmetry in — > K^gK^ into 
account as described in Refs. [3,[22l. Both the shapes and 
the normalizations of the misidentification backgrounds 
are fixed in the fit. The asymmetry and the sum of the 
D"*" and D~ yields are directly obtained from a simulta- 



To obtain Al* we first extract A^"""^ from 
a simultaneous fit with the same parameterizations 
for the signal except for the misidentification back- 
ground. The values of A^^^^ ^ ^ are evaluated in 
4x4x4x4x4 bins of the five-dimensional (5D) phase 
space {ptk-, cos6'^>i, pJp^+, cosS^^, cos^g^^). Each 
candidate is then weighted with a fac- 
tor of I =F^^(.^^ ^ ^ in the corresponding bin of the 5D 
phase space, where the phase space of the tt^ with lower 
Pt in — >■ K^TT^TT^ decay is used. After this weight- 
ing, the asymmetry in — s- K~Tr~^TT~^ decay sample 
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becomes , where tt"*" refers to the tt"*" with higher pT 
in the decay. The detector asynimetry, , is measured 
from simuhaneous fits to the weighted M^K^tt^tt^) dis- 
tributions in 10x10 bins of the 2D phase space {Px^+, 
cos0j^+) with the same parameterization used in D'^ 

K~'k^t:^ decays. Figure [2] shows the measured A^^ in 
bins of and cosSj^^ together with A^^^-^'^'^^'^^ 

for comparison. The average of A'^^ over phase space is 
(-h0.078 ± 0.040)%, where the error is statistical. 

Based on a recent study of the A-p [Si , we obtain the 
asymmetry in bins of Kg momentum in the lab. For the 
present analysis, Aj) is approximately 0.1% after inte- 
grating over the phase space of the two-body decay [l^l ■ 

The data samples shown in Fig. [T] are divided into 



qlab 

KgTT^ candidate is then weighted 



10x10x16 bins of the 3D phase space {Px^+, cos9 
p'^t^). Each D± 



with a factor of (1 =F A'^ )(1 T Ax)) in the 3D phase 
space. The weighted M{Kgn^) distributions in bins of 
cos 6'^^ are fitted simultaneously to obtain the corrected 
asymmetry. We fit the linear component in cos 9'^^ 
to determine ApB while the Acp component is uni- 
form in cos 

as a function of | cos 6'^^ \ . From a weighted aver- 
age over the |cos6'^'^^| bins, wc obtain A^p^^^^^ = 
(-0.363 ± 0.094)%, where the error is statistical. With 
out the A-p correction as in previous publications 
the value of Acp is (-0.462 ± 0.094)%. 
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FIG. 3: Measured Acp (top) and Afb (bottom) values as a 
function of | cosS^^^|. In the top plot, the dashed line is the 
mean value of Acp while the hatched band is the ilcrtotai 
interval, where atotai is the total uncertainty. 

The method is validated with fully simulated Monte 
Carlo events [i^ and the result is consistent with no in- 
put asymmetry. We also consider other sources of sys- 
tematic uncertainty. The dominant one in the Ac p mea- 



surement is the A^ determination, the uncertainty of 
which is mainly due to the statistical uncertainties in the 
D+ A'^7r+7r+ and K^n+n^ samples. These 

are found to be 0.040% and 0.048%, respectively, from 
a simplified simulation study. A possible Acp in the 
D'^ — > X^TT+TT^ final state is estimated with the relation, 
Acp — —y sin S sin 4>y/R [g^]- Using the 95% upper and 
lower limits on — Z)° mixing and CP violation pa- 
rameters [i^l, Acp in the D'^ — > K^tt^tt'^ final state is 
estimated to be less than 0.014% and this is included as 
one of systematic uncertainties in the A^ determination. 
By adding the contributions in quadrature, the system- 
atic uncertainty in the A^ determination is estimated 
to be 0.064%. We estimate 0.003% and 0.008% system- 
atic uncertainties due to the choice of the fitting method 
and that of the cos6'^^ binning, respectively. Finally, 
we add the systematic uncertainty in the Ap correction, 
which is 0.016% based on Ref. [18[. The quadratic sum 
of the above uncertainties, 0.067%, is taken as the total 
systematic uncertainty. 



We find A 



CP 



= (-0.363±0.094±0.067)%. This 
measurement supersedes our previous determination of 

[Zl. In Table H we compare all the available 



A 



CP 



measurements and give the new world average. 



TABLE I: Summary of A(-,p ^ measurements (where the 
first uncertainties are statistical and the second systematic), 
together with their average (where only the total uncertainty 
is quoted). 



Experiment 


+ (%) 


FOCUS [5] 


-1.6±1.5±0.9 


CLEO [6| 


-1.3±0.7±0.3 


BaBar [8] 


-0.44±0.13±0.10 


Belle (this measurement) 


-0.363±0.094±0.067 


New world average 


-0.41±0.09 



According to Grossman and Nir [l9|, we can estimate 
the experimentally measured CP asymmetry induced by 
SM K'^ — mixing, A^p, assuming negligible DCS de- 
cay D+ K°n+ in the final state D+ K°gTr+ . By 
multiplying Aj^p by the correction factor 1.022 ± 0.007 
due to the acceptance effects as a function of Kg decay 
time in our detector, we find the the measured asymme- 
try due to the neutral kaons to be (—0.339 ± 0.007)%. 

In summary, we report evidence for CP violation in 
the decay Z?+ ^s'^^ using a data sample correspond- 
ing to an integrated luminosity of 977 fb~^ collected with 
the Belle detector. The CP asymmetry in the decay is 
measured to be (—0.363 ± 0.094 ± 0.067)%, which repre- 
sents the first evidence for CP violation in charmed me- 
son decays from a single experiment and a single decay 
mode. After subtracting the contribution due to K^ — K^ 
mixing (A^p), the CP asymmetry due to the change 
of charm (App = A^p~^^°'^^) is consistent with zero. 
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A^C' ^ (-0.024 ±0.094 ±0.067)%. The measurement in 
the decay _D+ — > Kgir'^ is the most precise measurement 
of Ac p in charm decays to date and can be used to place 
stringent constraints on new physics models in the charm 
sector 
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